Current-voltage characteristics of graphene devices: interplay between Zener-Klein 

tunneling and defects 
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We report a theoretical/experimental study of current-voltage characteristics (I-V) of graphene 
devices near the Dirac point. The I-V can be described by a power law (I oc V a , with 1 < a < 1.5). 
The exponent is higher when the mobility is lower. This superlinear I-V is interpreted in terms of 
the interplay between Zener-Klein transport, that is tunneling between different energy bands, and 
defect scattering. Surprisingly, the Zener-Klein tunneling is made visible by the presence of defects. 
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Zener tunneling^ is a concept known since the 30's, 
which, in a solid, refers to the tunneling of carriers from 
one band to another through the forbidden energy gap 
(for example from the conduction to the valence band). 
This tunnel process is very intriguing in graphene be- 
cause the energy gap is suppressed to zero and because of 
the peculiar charge carriers behaving as Dirac fermions 
In particular, some of the carriers (those with the 
velocity parallel to the electric field) experience Zener 
tunneling without being backscattered — ~— , a behavior 
which is markedly different from the one in conventional 
semiconductors. The physics is the same as for relativis- 
tic electrons tunneling through a barrier, a phenomenon 
called Klein tunneling 7 and, for this reason, we will use 
the term Zener-Klein (ZK) tunneling. 

In view of the remarkable properties of ZK tunneling in 
graphene, it is understandable that an intensive endeavor 
was made to challenge it. So far, the effort was focused 
on graphene p-n junctions^—. In these devices, carri- 
ers tunnel through a sharp energy barrier induced with 
external local gate electrodes. Sophisticated nanofabri- 
cation techniques were employed to structure these local 
gates. For instance the insulator layer was very thir*^^ii^, 
the local gate was separated from the graphene by an air 
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or the local gate was extremely narrow^* 



gap 

Here, we argue that Zener-Klein tunneling can be ob- 
served in graphene with the most common device lay- 
out (undoped, four-point configuration, and without any 
local gates) by simply measuring the I-V at room tem- 
perature. First, we provide an analytical semi-classical 
expression for the I-Vs as a function of the doping. In 
graphene, the ZK current manifests itself with a super- 
linear current / oc V a , with a — 1.5. Then, we study the 
role of defects with the "exact" (non-perturbative) non- 
equilibrium Green-function approach finding the coun- 
terintuitive result that charged impurities enhance the 
visibility of the ZK current. Finally, we report measure- 
ments showing that the I-Vs at the Dirac point is indeed 
described by power laws, / oc V a , with a ranging from 
1 to 1.4. The exponent a is higher when the mobility is 
lower, consistently with our theoretical predictions. 

In graphene ZK tunneling leads to unusual I-Vs as 
compared to those of semiconductors/insulators. Let us 



consider transport through a piece of a material and ap- 
ply a voltage -V between the right (R) and left (L) sides. 
For a semiconductor with electronic gap E g , ZK tunnel- 
ing is possible only for eV > E g , where e > is the 
electron charge (Fig. [1]). On the contrary, in graphene 
(usually defined as a semi-metal) the gap is zero and, 
thus, ZK tunneling is possible for arbitrarily small V. 

More specifically, the two-dimensional electronic-band 
dispersion of graphene is a cone: e = ±hvp\/k± + k 2 , 
where k (fej_) is the wavevector-component parallel (per- 
pendicular) to the current flow. During ballistic trans- 
port (in absence of scattering) k± is conserved. For a 
fixed the bands are hyperbolae with gap A = 2hvpk± 
(Fig. [IJ. For any V, there are conducting channels 
for which the tunneling is possible (with k± such that 
A < V"). We will show that this results in a tunneling 
current I oc V 3 ^ 2 . By contrast, the ZK tunneling current 
in semiconductors vanishes exponentially at low V. 
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FIG. 1: (Color online) Bands of the L and R contacts in a 
semiconductor and in graphene. The arrows represent the 
possible occurrence of Zener-Klein tunneling. 

In graphene, within the Landauer approach, the cur- 
rent per unit of lateral length, J, is 
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T(e, V)de 
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where the factor 4 accounts for spin and valley degener- 
acy and the transmission T(e, kj_,V) is the probability 
that an electron (with energy e and perpendicular mo- 
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mentum fcjj is transmitted through the channel. We as- 
sume a uniform drop of the electrostatic potential along 
the current-flow direction, with constant electric field 
V/l, being I the distance between the contacts. 
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FIG. 2: (Color online) Electronic transport in graphene: the- 
ory, a: electronic bands of the left (L) and right (R) contacts. 
The hyperbolae are the bands corresponding to a finite k± 
(A = 2fiVFk± — 22 meV). In the L and R contacts the bands 
are filled up to cf and ef — eV, respectively, where £f > 
(ef < 0) corresponds to electron (hole) doping. b,c: electronic 
transmissions T and T, defined as in Eq.[T] for V — 0.1 V and 
I = 1 /im. In the gray zone (the region corresponding to the 
Zener-Klein tunneling) T and T are magnified, b: ballistic 
case calculated with NEGF or with the WKB approximation, 
c: NEGF results in the ballistic case (no defects) or in the 
presence of long-range (lr) or short-range (sr) defects. d,e: 
calculated current I per unit of lateral-length lw vs. V (the 
voltage applied between the electrodes) as a function of the 
gate voltage (V g ). V g goes from to 35 V with 5 V step, 
I = I ^m. In the ballistic case (d), lines are approximated 
semiclassical results (i.e. the analytical curves from noteQ3), 
points are from the "exact" NEGF simulations. 

The transmission can be calculated with the non equi- 
librium Green function (NEGF) method. To describe 
the purely ballistic case, we also use a semiclassical ap- 
proach, based on the Wentzel-Kramers-Brillouin (WKB) 
approximation. The transmission T WKB (e, k±, V) can 
be equal to 1, 0, or to T ZK = cxp[-nlA 2 /(4:hv F eV)p^. 
(see the example in Fig. [5Jib). We call non-tunneling 
current (Fig. [5Ji), the one associated with carriers that 
always remain in the same band n or tt* (T wkb = 1, 
light-shadowed (yelow) area in Fig. [5^,). We call "Zener- 
Klein" current, the one associated with carriers that tun- 
nel from the tt to the tt* band (T WKB = T ZK , dark- 
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FIG. 3: Electronic transport in graphene: measurements, a: 
measured current / per unit of lateral-length lw vs. V in 
a low-mobility sample, for different gate values (Vg=0,-20,- 
40 V). V g has been shifted to -14 V so that V g =0 corresponds 
to the Dirac point. The length between the voltage elec- 
trodes I — 1.1 jitm and lw = 1.1 /im. b: measured I-V in 
a high-mobility sample, at the Dirac point, before and after 
the introduction of defects through electronic bombardment. 
Z=2.2 lira and lw=550 nm. c: double- logarithmic scale plot 
of the I-V . d: exponent a as a function of mobility fi for 
different devices. I varies from 0.9 to 5.9 /im and lw from 70 
to 1500 nm. 



shadowed (cyan) area in Fig. [5^). From Fig. the 
rpWKB transmission is a good approximation to our most 
precise NEGF calculations^. 

In a graphene-based field-effect device, the density of 
the carriers n can be varied by changing the gate volt- 
age V g . n = VgCg/e, being C g the gate-channel capaci- 
tance. Fig.[5Ji reports the current-voltage (/-V) curves in 
the ballistic regime obtained with the semiclassical WKB 
approach (by letting T = T WKB in Eq. [TJ and with the 
"exact" NEGF method 15 , for various dopings (we use 
C g = 1.15 x 1CT 4 F/m 2 1S ). The two methods give al- 
most identical results (for the WKB case, we report in 
note [l^ an analytical expression for the I-V as a func- 
tion of ep). For zero-doping (V g = V) there is no 
contribution from the non-tunneling current, the current 
is entirely due to ZK tunneling, and the I-V curve is su- 
perlinear (I oc T/ 3 / 2 )!?. As soon as the system is doped 
(already for V g = 5 V) the ZK current is no more domi- 
nant (with respect to the non-tunneling current) and for 
small bias (V < 0.1 V) the I-V is ohmic (linear). 

Do we expect the superlinear ZK current to be visible 
in actual devices? At first sight the answer is no for two 
reasons. First, in actual devices, the carrier concentra- 
tion is never exactly zero. Indeed it has been observed2i 
that the presence of charge impurities induces a spatial 
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fluctuation of the Fermi level with respect to the Dirac 
point. As a result, it is difficult to achieve the experimen- 
tal condition where ZK tunneling is observable (V g =0 in 
Fig. [5Ji). Second, the scattering of the carriers with opti- 
cal phonons with energy fiw=0.15 eV causes the current 
to saturate when increasing V to high value d 18 ' 20 . This 
process occurs for eV > hwl/l e i (l e i is the carrier elas- 
tic scattering length, due to defects) and is, thus, par- 
ticularly relevant for high-quality high-mobility samples 
(with high l e i). This saturation of the non-tunneling cur- 
rent induces a sublinearity (/ oc V a , with a < 1) which 
tends to cancel the super linearity (a > 1 ) of the ZK cur- 
rent, further masking it (see [13 for further discussion). 

The situation is possibly changed by the presence of de- 
fects. Actual devices are characterized by defects which 
scatter electrons elastically (that is conserving the en- 
ergy)^. Elastic defects can be neutral point defects 
or charged Coulomb impurities^ outside the graphcnc 
plane (usually at a distance ~ 1—2 nm) — . Point de- 
fects affect the electrostatic potential seen by the carri- 
ers on a length scale smaller than the graphene unitary 
cell (short — range) and, thus, the carriers cannot be de- 
scribed in terms of electronic bands. On the other hand, 
charged impurities modify the potential uniformly on a 
length scale much longer than the unit cell (long — range) 
and the electronic bands are still a meaningful concept. 
The ZK current is expected to be more sensitive to short- 
range defects than to long-range ones. Indeed, the ZK 
current is determined by a transition between two bands 
whose relative energy is not affected by long-range de- 
fects. Moreover, long-range defects are expected to di- 
minish the non-tunneling current. Overall, one could 
wonder whether the presence of long-range defects can 
be used to suppress the non-tunneling current and, thus, 
to make visible the ZK one. 

To verify this hypothesis, we simulate disordered 
graphene within NEGF by considering both long- and 
short-range elastic defects^. We remark that the NEGF 
approach provides an exact (non-perturbative) atomistic 
treatment of disorder. Defects are simulated by changing 
randomly the on-site potential by Vd — 0.1 eV. This Vd is 
a realistic choice since it provides a low-bias conductivity 
in reasonable agreement with measurements^. 

From NEGF simulations, the presence of long- 
range defects diminishes the non-tunneling transmission 
(Fig. [2b) but, in general does not reduce the ZK one. 
For V—0.1 V, long-range defects even increase the ZK 
transmission (Fig. 120 . We checked that short-range de- 
fects diminish, as expected, both the non-tunneling and 
the ZK transmission, with respect to the ballistic case 
(Fig. [2b)- To see whether the relative increase of the ZK 
transmission with respect to the non-tunneling one can 
lead to measurable effects, in Fig. [2b we show the theo- 
retical I-V curves in the presence of long-range defects. 
The superlinear behavior (the signature of the ZK cur- 
rent) is still visible at V g = (/ oc V a , with a = 1.4 in 
Fig. [2b) and is also visible at finite . 

We now turn our attention to measurements, carried 



out on single- layer graphene devices^. Different devices 
were fabricated in a four-point configuration and have dif- 
ferent mobilities [i ranging from 80 to 20000 cm 2 V _1 s _1 
(low mobility corresponds to a higher density of de- 
fects)^. Fig. [3k shows a typical set of I-V characteristics 
for different V g applied on the backgate for a sample with 
a relatively modest mobility (/x = 1700 cm 2 V — 1 s — •*•). The 
I-V is superlinear at the V g of the Dirac point, consistent 
with the above prediction of ZK tunneling. The super- 
linearity is better seen in a double-logarithmic scale plot 
(Fig. [3b) where the I-V is reasonably well described by a 
power law I oc V a with a = 1.3. Both a and the current 
values are in a remarkable qualitative agreement with cal- 
culations given the simplicity of the model as can be seen 
by comparing Fig. [2]3 (Z=l/im) and Fig. [3Ji (Z=l.l/xm) 
for small V g . More elaborated models (e.g. with a more 
realistic description of impurities and including electron- 
phonon scattering) are required to reach a quantitative 
agreement between theory and measurements. 

We observe that the superlinearity vanishes for devices 
with high /i. Fig. [5Ji shows a (extracted at the Dirac 
point) as a function of the mobility /i of 22 different de- 
vices. Indeed, as the mobility increases, a tends to 1 (cor- 
responding to linear I-V). In an additional experiment, 
we introduced defects in a high-mobility graphene device 
by bombarding it with 10 keV electrons. From Fig. [3fc>, 
before bombardment the mobility n =7000 cm 2 V _:l s — 1 
and the I-V is linear with a — 1.0. After bombard- 
ment \x drops to 260 cm 2 V _1 s _1 and the I-V becomes 
superlinear (a = 1.2). These observations are consistent 
with the above discussion that in high-mobility samples 
the ZK superlinearity is masked by the non-tunneling 
current. Namely, the reduction of disorder increases the 
contribution of the non-tunneling current with respect to 
the ZK one and, also, favors the non-tunneling current 
saturation (due to scattering with optical phonon o 18 ' 20 ) . 

We now discuss other mechanisms that could lead to 
superlinear /-Vs. It could be related to the physics oc- 
curring at tunnel barriers (such as the Luttinger liquid- 
like behavior in nanotubes or the breakdown of insu- 
lating barriers). However, measurements are done on 
high-quality devices in a four-point configuration, which 
makes the presence of tunnel barriers unlikely. Superlin- 
ear I-Vs could also be attributed to quantum effects, such 
as weak localization or electron-electron interaction, but 
these effects should be negligible since the applied current 
is large, heating the graphene layer to several hundreds 
of Celsius 2 ^. Overall, Zener-Klein tunneling remains the 
most plausible mechanism to explain our measurements. 

We finally stress that previous observations of Klein 
tunneling 10 ' 12 i 13 in graphene were done using very dif- 
ferent device setups. In [l0lfl2Tfl"3l the carriers tunnel from 
conduction to valence bands in a p-n junction. In these 
nanostructured devices, the ZK tunneling is observed 
thanks to a configuration which allows to eliminate the 
non-tunneling current and thanks to the intense electric 
field at thep-n junction (~ 10~ 3 eV/A, see suppl. info 
of [l3l and Hf). In our devices, which are not p-n junc- 
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tions, the non-tunneling (intraband) current is present 
(this current can mask the ZK tunneling one) and the 
electric field (~10~ 5 eV/A) is substantially weaker. De- 
spite these unfavorable conditions, it is possible to probe 
the Zener-Klein effect. 

Concluding, measurements and calculations show, con- 
sistently, that the I-Vs of graphene devices become su- 
perlinear in the presence of disorder (in low-mobility sam- 
ples). The superlinearity is attributed to Zener-Klein 
tunneling (tunneling between different energy bands, 



from 7r to 7T*). In high- mobility (high-quality) graphene 
samples, the superlinearity is masked by the contribu- 
tion of the non-tunneling current (due to carriers al- 
ways remaining the same band). In low- mobility sam- 
ples, the Zener-Klein tunneling current is visible because 
the higher density of defects decreases (filters) the non- 
tunneling current. 
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